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Abstract: To accelerate the fractal decoding process, a minimum iterated function system based fast 
fractal decoding method was proposed in this study. In fractal encoding process, we found that there exists 
a minimum domain block set (MDBS) which can provide the best-matched domain blocks for all range 
blocks, and then the domain blocks of MDBS, the range blocks within MDBS, and the associated mapping 
operations between them constitute the minimum iteration function system (MIFS). In decoding process, 
MIFS was first recovered in the first iteration. Then, in each of the second to penultimate iterations, only 
the range blocks within MDBS are reconstructed, and the computations of reconstructing the remaining 
range blocks can be saved. Finally, all range blocks are reconstructed to obtain the decoded image in the 
last iteration. Four fractal encoding methods were adopted to assess the performance of the proposed 
method. Experimental results show that the proposed method can complete the decoding process with 
fewer computations. 
Keywords: Fractal image decoding; Minimum iterated function system; Minimum domain block set; 
Redundant domain block set; fixed point 
1. Introduction 

The kernel idea of fractal image coding was first proposed by Barnsley, and then Jacquin 
proposed the first practical fractal coding algorithm which was considered as a promising 
image compression method!'?!, Although suffering from high computational complexity in 
encoding process, fractal image coding has the characteristics of novel idea, potential high 
compression ratio, fast decoding, and resolution independence. Thus, on one hand, to 
accelerate the encoding process, researchers worldwide proposed all kinds of fast fractal 
encoding methods which can be mainly divided into two categories: 1) Local block matching 
based fast fractal image coding]. By converting exhaustive block matching into local block 


matching in domain block pool, this type of methods tried to accelerate the encoding process 


while maintaining the decoded image quality as much as possible. 2) No-search fractal image 
coding!**!, This type of methods appointed the best-matched domain block directly without 
block matching operations and can realize real-time encoding at the expense of poor decoded 
image quality. On the other hand, after many years’ research, fractal image coding has been 


gradually applied to many other image processing fields, such as image denoising!'*'®], image 


[17-21] [22,23] 24-26] [27-29 


magnification , image hashing , image retrieval!**+76l, watermarking!*’"!, and human 
pose estimation?°3!], 

In our study, a novel method was proposed to accelerate the fractal decoding process. 
Fractal encoding process generally consists of a series of mapping operations from domain 
blocks to range blocks. First, we found that the block matching operations can be completed 
only with partial domain blocks which are defined as the minimum domain block set 
(MDBS). Then, the remaining domain blocks outside MDBS are defined as the redundant 
domain block set (RDBS). Further, the domain blocks of MDBS, the range blocks within 
MDBS, and the associated mapping operations are defined as the minimum iterated function 
system (MIFS). From the above definitions, we know that by excluding the domain blocks of 
RDBS, MIFS contains necessary mappings operations from MDBS to all range blocks and 
plays a crucial role in encoding process. In decoding process, in first iteration, all range 
blocks are reconstructed. By distinguishing MDBS from RDBS, MIFS can be recovered. In 
each of the second to penultimate iterations, only the range blocks within MDBS are 
reconstructed by MIFS. In the last iteration, the local image within MDBS converges, and 
then the range blocks within RDBS are reconstructed as well. In summary, in the first and last 
iterations, we perform the same operations with the previous method. But for each of the 
second to penultimate iterations, because the domain blocks within RDBS provide no 
mapping operation for any range block, there is no need to reconstruct the range blocks 
within RDBS, and the associated computations can be saved. Three state of the art methods 
and Jacquin’s method are adopted to assess the performance of the proposed method. 
Experimental results show that compared with the conventional method, the proposed method 
can effectively reduce the number of computations in decoding process. 


2. Conventional fractal image coding 


Conventional fractal image encoding proposed by Jacquin can be described in detail as!!: 
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The M x N input image is first divided into nonoverlapping B x B range blocks, 


R,,i=1,2,3,...,NumR , where NumR denotes the total number of range blocks. Then, the 


domain blocks, D;,/=1,2,3,..., NumD , can be obtained by sliding a 2B x 2B window over 


the input image with a fixed sliding step 6=2B, where NumD denotes the total number of 
domain blocks. Further, all domain blocks are contracted to the size of B x B and constitute a 
domain block pool which is then extended with eight isometric transformations. Finally, for 
arbitrary range block Ri, the best-matched domain block D; can be found by minimizing the 


following function: 


CE(R,)=argmin | 
D; 


R, — Q; (D, )|- as le - fay, [r p, Ji AL | 


i=1,2,3,..., NumR, j =1,2,3,..., NumD 


(1) 


where 9%; (-) denotes the mapping operation from PD; to R,, n (-) and y;(-) denote the 
contracting and isometric transformations, respectively, and @; and Ø; denote the scaling 
and offset coefficients of the affine transformation, respectively. J denotes the B x B block 


whose components are all ones. 


(d) Third iteration (e) Fourth iteration (f) Fifth iteration 
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Fig. 1 Initial image and the first five iteration images in decoding process. 

In decoding process, arbitrary M x N image can be selected as the initial image. In each of 
the following iterations, based on the same operations in encoding process, all range blocks 
are reconstructed by their respective best-matched domain blocks. After about ten iterations, 
the decoding process converges to the decoded image. When the Bridge and Lake images are 
selected as the input image in encoding process and initial image in decoding process, 
respectively, Figure 1 (a) illustrates the initial image, and the remaining five images in Fig.1 
(b)-(f) illustrate the first five iteration images in decoding process. We can see that the initial 
image can gradually converge to the input image. 


3. Accelerating fractal decoding process by minimum iterated function system 


(c) Mappings of D, (d) No Mapping provided by D, 


Fig. 2 Illustration of mapping operations from domain blocks to range blocks. 


In encoding process, by a series of contracting, isometric, and affine transformations, the 


range blocks can be approximated by their respective best-matched domain blocks as: 


Ü e (wÜ fez [n (0) a4] 


i=1,2,3,..., NumR, j =1,2,3,..., NumD' 


NumR 
I= R x 
U : (2) 


i= 


where D; denotes the best-matched domain block of R,, NumD' denotes total number of 
the best-matched domain blocks we actually need and satisfies NumD’ < NumD . Although 


all domain blocks, Dj, j =1,2,3,...NumD , are the candidates of the best-matched domain 


block of R, , only partial domain blocks, Dj, j =1,2,3,...,NumD’ , are the ones we actually 


need. The remaining domain blocks cannot provide mapping operation for any range block 


and will be not used at all. For example, in Fig.2 (a), we have three domain blocks, 


D,,j =1,2,3 , and each domain block can be uniformly divided into four range blocks, 


R,,i=1,2,3,4. The arrow “—” represents the mapping operation 9, (-) from one domain 
block to another range block. Figure 2 (b), (c), and (d) illustrate the mapping operations that 


are provided by domain blocks D,, D,, and D;, respectively. We can see that in Fig.2 (b) 
and (c), both D, and D, can provide mapping operations for range blocks. But in Fig.2 (d), 


D, provides no mapping operation for any range block. Thus, we can divide the domain 


blocks into two categories: Minimum domain block set (MDBS) and redundant domain block 
set (RDBS): 
Definition 1: The domain blocks that can provide the mapping operations for range blocks 


more than or equal to once constitute the minimum domain block set (MDBS): 


MDBS =|_)D’ s.t. Times(Dj, R,) > 1,/=1,2,3,..., NumR, j=1,2,3,....NumD' 6) 
J 


where Times(A, B) denotes the total number of mapping operations from A to B. From 
Definition 1, we know that the domain blocks of MDBS can provide mapping operations for 
all range blocks (including the range blocks within MDBS), In Fig.2, both D, and D, 


represent the domain blocks of MDBS, and then we have the definition of redundant domain 


block set (RDBS) as: 
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Definition 2: The domain blocks that provide no mapping operation for any range block 
constitute the redundant domain block set (RDBS): 

RDBS = UP; s.t. Times (D", R = 0, i=1,2,3,..., NumR, j7=1,2,3,. .., NumD' (4) 
where NumD" denote the total number of the domain blocks of RDBS and satisfies 
NumD"=NumD —- NumD'’ . In Fig.2, D, represents the domain block of RDBS. For Lake 
and Bridge images in Fig.3 (a) and (b), the domain blocks within the white and red boxes 


represent the domain blocks of MDBS and RDBS, respectively. 
se oH 5 Geocmasee neu a PREECE EEFE i 
im fe +t 


(a)Lake (b) Bridge 
Fig.3 Illustration of MDBS and RDBS. 
Further, we have the definition of minimum iterated function system (MIFS) as: 
Definition 3: The domain blocks of MDBS, the range blocks within MDBS, and the 
associated mapping operations between them are defined as the minimum iterated function 
system (MIFS): 


MIFS:|(R},D‘),g, | i=1,2,3,...,NumR’, j =1,2,3,..., NumD', D; €eMDBS (5) 


where R/,i=1,2,3,...,.NumR’ denote the range blocks within MDBS. NumR’ denotes the 


total number of the range blocks within MDBS and satisfies NumR’=NumD'’x4 . From 
Definition 1, 2, and 3, we know that because MIFS can provide mapping operations from the 
domain blocks of MDBS to the range blocks within MDBS, MIFS can complete the mapping 
operations by itself independently. Moreover, because the domain blocks of MDBS within 
MIFS can provide mapping operations for the range blocks within RDBS as well, MIFS can 
be used to reconstruct the whole decoded image in decoding process. Figure 4 illustrates the 


MIFS contained in Fig.2 (a). 


Fig.4 MIFS contained in Fig.2 


In decoding process, the proposed decoding method can be designed as: In first iteration, 
based on the mapping operations in encoding process, we reconstruct all range blocks, and 
MIFS can be recovered by distinguishing MDBS from RDBS. Then, in each of the second to 
penultimate iterations, based on MIFS, only the range blocks within MDBS are reconstructed 
by the domain blocks of MDBS. On one hand, for the range blocks within MDBS, MIFS 
maintains the same reconstruction operations as the previous method. Thus, we have the 
same local iteration image within MDBS as before which can be used to judge whether 
convergence has been achieved, and the same number of iterations can be also maintained. 
On the other hand, based on Definition 2, because the domain blocks of RDBS do not take 
part in reconstructing any range block, there is no need to reconstruct the range blocks within 
RDBS, and the associated computations can be saved. Finally, in the last iteration, after the 
range blocks within MDBS are reconstructed, the range blocks within RDBS are 
reconstructed as well, and the whole decoded image can be obtained. Figure 5 illustrates the 
process of the proposed decoding method, and we can see that the proposed method maintain 
the same number of iterations as the previous method. In the first and last iterations, the 
unshaded regions represent all range blocks within MDBSs and RDBSs we need to 
reconstruct, and the proposed method does the same operations with the previous method. 
For each of the second to penultimate iterations, the unshaded regions still represent the range 
blocks within MDBSs we need to reconstruct. However, the shadow regions represent the 
range blocks within RDBSs whose reconstruction operations can be saved, and thus the 


decoding process can be accelerated. 
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Fig.5 Flow chart of the proposed decoding process. 


4. Experiments 

In this section, two 256 x 256 images, Lake and Bridge, are selected as the test images. The 
size of range blocks is set at 4 x 4, and the sliding step is 8. The scaling and offset 
coefficients, s and o, are quantized with 5 and 7 bits, respectively. Jacquin’s method and three 
state of the art methods, Jacquin’s’!, Chauraisa’s!?!, Zheng’s!“!, and Gupta’s methods!’!, are 
adopted to assess the performance of the proposed method. The root of mean square error 
(RMSE) is adopted to measure the deviation between the kth iteration and decoded images as: 


LESS. p (./) 


j=l i=l 


RMSE“) = k=1,2,...,N (6) 


Hxw 


where f and fPeeded denote the Ath iteration and decoded images, respectively, H and W are 
the height and width of the images, respectively, and N denotes the total number of iterations 
we need to complete the decoding process. In the experiment, the decoded image can be 
obtained by encoding and decoding the input image in advance, and the detailed experimental 
procedures are listed as: 

Step 1: Select one fractal encoding method and encode the input image. 

Step 2: In the first iteration of decoding process, select a 256 x 256 blank image as the 
initial image, reconstruct all range blocks, and all domain blocks can be divided into MDBS 
and RDBS. Let k=1, calculate and record RMSE within MDBS. 

Step 3: Reconstruct the range blocks within MDBS. Let k=k+1, calculate and record 
RMSE within MDBS. If the local image within MDBS satisfies the convergence 
requirement of RMSE") -RMSE“”|<1 , go to Step 4. If not, turn back to Step 3, and 
perform another iteration. 

Step 4: Reconstruct the range blocks within RDBS, and the final decoded image can be 


obtained. 


In Tables 1, 2, 3, and 4, “A” represents achieving convergence in the Mth iteration. For 
Lake image in Tab.1, we totally have NumD=1024 domain blocks, and MDBS and RDBS 
contain NumD’=922 and NumD"=102 domain blocks which comprise 


NumD'/NumD x100%=90.04% and NumD"/NumD x100%=9.96% of the input image, 


respectively. In each of the second to (N-1)th iterations, we only reconstruct the range blocks 
within MDBS. On one hand, the range blocks within MDBS are reconstructed by the domain 
blocks of MDBS, which implies that the same reconstructing operations are performed within 
MDBS as Jacquin’s method. Thus, the proposed method can keep the same RMSEs within 
MDBS as Jacquin’s method, and the same number of iterations can be also maintained as 


Jacquin’s method. Then, when the deviation between the current and previous local iteration 
images within MDBS achieves the requirement of RMSE“) -RMSE“™”|<1 , the decoding 


process can be considered as convergence. On the other hand, because the proposed method 
only reconstruct the range blocks within MDBS rather than all range blocks in Jacquin’s 
method, 9.96% of total computations can be saved regarding Jacquin’s method. Finally, all 
the range blocks within both MDBS and RDBS are reconstructed in the last iteration. For the 
whole decoding process, the proposed method does the same with Jacquin’s method in the 
first and last iterations. For each of the second to (N-1)th iterations, only the range blocks 
within MDBS are reconstructed. In each iteration, if the percentage of the computations 
required (PCR) for reconstructing all range blocks is considered as 100%, the PCR for the 
proposed method in the whole decoding process can be calculated as: 


_ (N-2)x(NumD'/NumD )x 100% + 2 x 100% 
p N x100% 


PCR x100% (7) 


where the numerator and denominator represent the computations required in the proposed 


and previous methods, respectively. Based on Eq. (7), we know that for Lake image in Table 
1, we have NumD'/NumD x100%=90.04% , and N=11. Thus, PCR=91.85%, and 8.15% of 
total computations can be saved regarding Jacquin’s method. Similarly, for Bridge image, 
MDBS and RDBS contain NumD'=893 and NumD"=131 domain blocks, respectively, 
and N=11. By Eq. (7), we have PCR=89.54%, and 10.46% of total computations can be 


saved. 


Tables 2, 3, and 4 list the experimental results for Chaurasia’s, Zheng’s, and Gupta’s 


methods, respectively. In Table 2, for Lake and Bridge images, 81.35% and 79.98% of total 


computations are required in each of the second to ninth iterations, and PCRs are 84.74% and 


83.62%, respectively. In Table 3, for the two test images, 90.14% and 87.21% of the total 


computations are required in each of the second to tenth iterations, and PCRs are 91.93% and 


89.54%, respectively. In Table 4, for the two test images, 90.72% and 88.09% of total 


computations are required in each of the second to tenth iterations, and PCRs are 92.41% and 


90.26%, respectively. In summary, compared with the previous methods, the proposed 


method can effectively reduce the number of computations in decoding process. 


Table 1 Performance comparison between Jacquin’s and proposed methods in decoding process!” 


Iterations (k) PCR 
Image Methods Performance 
1 2 3 4 5 6 7 8 9 10 11 (%) 
j — RMSE(MDBS) 97.18 67.50 46.56 31.01 20.33 13.21 852 545 344 2.14 1.29 x 
acquin’s 
=) i PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 
ake 
E J RMSE(MDBS) 97.18 67.50 46.56 31.01 20.33 13.21 8.52 545 344 2.14 1.294 x 
ropose 
j PCR(%) 100 90.04 90.04 90.04 90.04 90.04 90.04 90.04 90.04 90.04 100 91.85 
J iG RMSE(MDBS) 81.66 55.01 37.56 24.95 16.62 11.02 7.26 4.74 3.06 1.94 1.194 x 
acquin’s 
ea i PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 
ridge 
j p d RMSE(MDBS) 81.66 55.01 37.56 24.95 16.62 11.02 7.26 4.74 3.06 1.94 1.194 x 
ropose 
E PCR(%) 100 87.21 87.21 87.21 87.21 87.21 87.21 87.21 87.21 87.21 100 89.54 
Table 2 Performance comparison between Chaurasia’s method and the proposed methodľ!. 
Iterations (k) PCR 
Image Methods Performance 
1 2 3 4 5 6 7 8 9 10 (%) 
Chaurasi RMSE(MDBS) 85.59 54.51 34.66 2146 13.08 7.87 4.70 2.79 1.64 0.954 x 
a a’s PCR(%) 100 100 100 100 100 100 100 100 100 100 100 
e 
7 d RMSE(MDBS) 85.59 54.51 34.66 2146 13.08 7.87 4.70 2.79 1.64 0.954 x 
ropose 
p PCR(%) 100 81.35 81.35 81.35 81.35 81.35 81.35 81.35 81.35 100 84.74 
Chaurasi RMSE(MDBS) 78.82 50.92 32.95 20.22 12.15 125 4.28 2.51 1.45 0.83 A x 
Brid a’s PCR(%) 100 100 100 100 100 100 100 100 100 100 100 
ridge 
B p d RMSE(MDBS) 78.82 50.92 32.95 20.22 12.15 7.25 4.28 2.51 1.45 0.83 A x 
ropose 
p PCR(%) 100 79.98 79.98 79.98 79.98 79.98 79.98 79.98 79.98 100 = 83.62 
Table 3 Performance comparison between Zheng’s method and the proposed method", 
Iterations (k) PCR 
Image Methods Performance 
1 2 3 4 5 6 7 8 9 10 11 (%) 
RMSE(MDBS) 97.38 67.76 46.81 31.24 20.51 13.35 8.62 552 3.49 2.17 131A x 
Lake Zheng’s 
PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 


Faos RMSE(MDBS) 97.38 67.76 46.81 31.24 20.51 13.35 8.62 5,92 3.49 2.17 1.314 x 
PCR(%) 100 90.14 90.14 90.14 90.14 90.14 90.14 90.14 90.14 90.14 100 91.93 
Ziengs RMSE(MDBS) 81.75 54.94 37.24 24.52 16.21 10.70 7.03 4.57 2.94 1.85 1.144 x 
iiiz PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 
Pora RMSE(MDBS) 81.75 54.94 37.24 24.52 16.21 10.70 7.03 4.57 2.94 1.85 1.144 x 
PCR(%) 100 87.21 87.21 87.21 87.21 87.21 87.21 87.21 87.21 87.21 100 89.54 
Table 4 Performance comparison between Gupta’s method and the proposed method"), 
Iterations (k) PCR 
Image Methods Performance 
1 2 3 4 5 6 7 8 9 10 11 (%) 
Canis RMSE(MDBS) 98.55 69.49 48.59 32.63 21.60 14.16 9.22 5.95 3.80 2.38 1.454 x 
Lake PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 
goes RMSE(MDBS) 98.55 69.49 48.59 32.63 21.60 14.16 9.22 5.95 3.80 2.38 1.454 x 
PCR(%) 100 90.72 90.72 90.72 90.72 90.72 90.72 90.72 90.72 90.72 100 92.41 
, RMSE(MDBS) 81.10 54.62 37.22 24.78 16.61 11.15 7.42 4.89 3.17 2.02 1.244 x 
Ree owe PCR(%) 100 100 100 100 100 100 100 100 100 100 100 100 
AE RMSE(MDBS) 81.10 54.62 37.22 24.78 16.61 11.15 7.42 4.89 3.17 2.02 1.244 x 
PCR(%) 100 88.09 88.09 88.09 88.09 88.09 88.09 88.09 88.09 88.09 100 90.26 


5. Conclusion 


A novel fast fractal decoding method is proposed in this paper. In encoding process, the 
definitions of MDBS, RDBS, and MIFS are introduced. In decoding process, through the 
same number of iterations as before, MIFS can complete the decoding process and achieve 
convergence within MDBS independently, and thus the computations of reconstructing the 
range blocks within RDBS from the second to penultimate iterations can be saved. Moreover, 
the reconstruction of the range blocks within RDBS in the first and last iterations should be 
performed to recover MIFS and obtain the whole decoded image, respectively. Finally, 
Jacquin’s and three state of the art methods are adopted to verify the effectiveness of the 
proposed method. 
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